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We present a global approach of non-dissipative physics. Based on symplectic mechanics this 
technique allows us to obtain the solution of a very large class of problems in terms of a Taylor 
expand. We apply this method to the problem of gravitational instability and we obtain a general 
expression of the gravitational potential, solution of the Vlasov-Poisson system, as a function of 
time in the context of Newtonian dust cosmology. 



I. INTRODUCTION 



The Hamiltonian formalism has been proved since a long time to be very general and efficient in classical dynamics. 
In the recent years, Hamiltonian techniques have also been applied efficiently to continuous systems. Here we give 
a general presentation of Hamiltonian formalism (reflecting the underlying symplectic structure) applied to systems 
with a finite or infinite number of degrees of freedom. The latter case means systems described by a distribution 
function in phase space (like electro-magnetic or gravitational plasmas). Exploring evolution in time, i.e., motion, 
we are able to derive, for such systems, the analytical expression giving the evolution of any physical quantity as a 
function of initial conditions only. This solution is given in terms of a Taylor development, with a rule to calculate 
the coefficients. 

Collisionless dynamics describes various astrophysical systems, typically self gravitating collapsing clouds of stars or 
dust particles. It may be applied, for instance, to study the stability and evolution of globular clusters, the formation 
of various cosmic structures by gravitational instability. Here we apply our Hamiltonian formalism to the study of 
^ gravitational instability, i.e., the evolution of a cosmic fluid under the influence of its own weight. Although our 
O calculations are performed in the Newtonian approximation, they apply perfectly to the evolution of perturbations 
j / - to the dust Friedmann - Lemaitre universe, assumed to describe the first stages of galaxies or large scale structure 
formation. For such systems we are able to derive the expression of any physical quantity as a development with 
respect to time, depending on initial conditions only. We illustrate by developing the value of the (self-consistent) 
gravitational potential at any point and at any time, as a function of initial conditions only. 

The paper is organized as follow : a first part is devoted to the presentation of the Hamiltonian formalism. 
Considering first systems with a finite number of degrees of freedom, we derive the evolution equation for any quantity, 
t j and explicit its solution. In order to show the mechanism of our formalism, we illustrate it on the very simple example 
of the harmonic oscillator. Then we generalize the formalism to continuous systems, by using functionals and functional 
i-H derivatives, which allows us in this case also to derive the equation evolution for any quantity, as well as its solution 
under the form of a development. This allows in principle to calculate any physical quantity as a function of initial 
conditions. 

H The second part specializes the problem to the gravitational instability. We first recall the dynamical equations de- 
f> i scribing it, and their solution in the cosmological context (Newtonian approximation to describe the dust Friedmann- 
Lemaitre models and their perturbations). Then we explicitly write the Poisson-Vlasov equations which describe the 
evolution of the system in the phase space, and its link with the usual hydrodynamical approach. Applying the results 
of the first part, we express the time-derivatives of physical quantities : as expected we recover the Liouville equation 
for the distribution function itself, and the derivative of the potential. We then explicit the Taylor development for the 
gravitational potential, and provide a recurrence formula which allows its explicit calculation from initial conditions. 



^1" 

(N 



i 

O 

5— ( 
^— » 



II. THE SYMPLECTIC APPROACH OF NON DISSIPATIVE PHYSICS 



This section presents a new and very general formulation of the non dissipative physical equations, applicable to 
many kinds of different problems. The unification proposed here is a generalization of several different works ( [1], 
[6]). We split the framework of non dissipative physic into two parts : systems with a finite number N of degrees of 
freedom on one hand, and we detail the case N = 1; on the other hand, systems with an infinite number of degrees 
of freedom, i.e., statistical physics. We will show that the equations for these two domains are strictly equivalent and 
derive from the same least action principle. We will propose moreover, the solution of this equation. 
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A. Systems with N degrees of freedom 



We will consider the 3-dimensional case, corresponding to one particle in 3-dimensional space. The variables are the 
spatial position q, and its conjugated quantity, the linear momentum p, both 3-vectors. The generalization to more 
dimensions is straightforward, the variables becoming any arbitrary collections of conjugated vectors (for systems of 
N particles for instance), or tensors. 

The foundation of our analysis relies on the possibility to associate a generator, a Hamiltonian-like function g(q, p), 
to any kind of transformation undergone by the system. In the following we will be mainly interested by one peculiar 
kind of transformation, i.e., motion. But transformations can also correspond to the effect of a class of perturbations 
on a given equilibrium state, or represent a geometrical change applied to the system (Rotation, Translation, etc.). 
In this sense, the formalism introduced here allows more than the simple study of motion : the search for invariances, 
the stability studies and, more generally, the systemization of the study of dynamical systems. 

To the transformation considered is associated a parameter A which allows to follow the transformation. When motion 
is the transformation considered, g is the usual Hamiltonian representing the energy E of the system, and A identifies 
with the time t (to study instabilities, g will represent the characteristics of the instability, and A the intensity of the 
displacement). Then one can always define the meta-action as the function 

*(p,q) = J ( p -^ _fi, ) rfA ' W 

where q and p are functions of A. In the case of motion, this meta-action reduces to the usual action, solution of 
the classical Hamilton-Jacobi equation (see [1] for example). In this peculiar case the total energy E of the system is 
conserved through the t— evolution. Similarly, our formulation requires that the generator g does not depend explicitly 
on A (although it depends through q and p) and is conserved through the transformation. 

Before turning to the key propositions, we recall some a standard definition and notations. The (usual) Poisson 
brackets for 2 functions x(p,q) and y(p, q) is defined as 

A / dx dy dx dy\ 

L 3 [x,y] := Jp 3 [x,y]dX, (3) 

which has all properties of a standard Lie bracket. 

We now state the fundamental proposition, whose demonstration follows from the calculations below : the evolution 
of p and q, as a function of A, is given by the least action principle 

L 3 [k,s] = [ P 3 [k,s] d\ = 0. (4) 



and we also define the new bracket 



for any function fc(q, p). This relation was initially presented in the special case of motion by [6] and [7]. 

Since (4) is true for any function k, it implies in particular Vi , Jj- = and Vi , Jj- = 0, corresponding to the two 
evolution equations for q and p : 



Vi ^- = ^ #- / (p.^-g)dX = 



• Firstly, 

dpi dpi J \ d\ 

dqi _ 9g_\ dx _ q 

d\ dpi d\ dpi J 

Since q and p are explicitly independent, and q.p vanishes on the surface of the system 1 , integration by parts 
gives 



J This assumption makes no problems in this finite dimensional case, more details are required for continuous systems, see 
section C and appendix A 



2 



<ki _ dg_\ dX _ 
dX dpi J 

which implies 

<hi_ _ dg_ _ d<udg_ _ dqi_dg_ rfq _ p 

dX dpi dqi dpi dpi dqi d\ ' ' 

which is the evolution equation for q induced by the transformation g. 
• Similarly, 

c% c% J \ dX J 



_d_<H _dy_\ dX _ 

' % dX dqi 



And, after a new integration by parts, we get 



dpi _ _dg_ _ dpi_dg_ _ dpi_dg_ dp _ ^ 

dX dqi dqi dpi dpi dqi dX ' ' 

the evolution equation for p induced by the transformation g. 
In the special case of motion, where g = E and A = t, (5) and (6) are Hamilton's equations. 



B. Evolution equation for any function 

More generally, the same technique allows to write the evolution equation of any function fc(q, p). Starting from 
the usual differentiation formula 

dk f dk dqi dk dpi \ 

~ ^ \dq~i~dX + dpi ~d~X ) ^ 

using (5) and (6) some straightforward algebra gives the evolution equation for k induced by the transformation g 

£ = ftlM- W 

This is a linear and first order differential equation. Its solution is known to be (when g is explicitly independent 
of A) 

A 2 

*(q, p) = *(q, p) k -A P 3 [g, k] \ Xo +y P 3 [ g , P 3 [ g , k}} \ Xo + ■ ■ ■ , (9) 

which is usually abbreviated under the notation 

fc(q,P) = e AP3[ ' 3] A(q,p) U„ • (10) 

This formulation is similar to the usual Hamiltonian formalism, with the main additional result here that explicit 
solution for the motion (or, more generally, the transformation) is given as a function of time (or, more generally A). 
This offers several advantages in comparison with the classic formulation. 

Firstly, we can expand any function in terms of any parameter. Although the case of motion is specially interesting, 
we have more general results. For example, taking for g the generator of a perturbation, and for k the energy of the 
system, we get a general energy variational principle for every system. In the context of stellar dynamics, [2], [8], [3] 
and [4] have produced new stability results. Finally, one can use this method to produce an integral of the evolution 
equation. A very important improvement lies in the fact that this formalism can be extended to systems with an 
infinite number of degrees of freedom. Before turning to it, we use a very simple application to illustrate how it works. 
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1. Application to the harmonic oscillator 



To illustrate some of the mechanism of the technique, we will shortly analyze here the trivial case of the motion of 
an one dimensional harmonic oscillator, with A = t and g = E = ^ + ^ui 2 q 2 . We use (9) with k =| q |:= q, and the 
initial conditions q(t = 0) = q and p(t = 0) = p , we have 

q(t) == q(t ) - tP 3 [E, q] \ t=to + ^P 3 [E, P 3 [E , q]] \ t=to -J^E, P 3 [E , P 3 [E , q]]] \ t=to +■■■■ (11) 

The Poisson brackets are easy to compute : 

P 3 [E q] = = _P_ 

' dq dp dp dq m 

P 3 [E,P 3 [E,q]] = P 3 [E,-^} = - q — 

m m 



P 3 [E,P 3 [E,P 3 [E,q}}} = P 3 [E,-^} = ^ 

and so on. We can obtain a recurrence rule. Injecting these results in (11) one can have 

p 1 q u> 2 2 t 3 p u> 2 t 4 q u> 4 t 5 p u> 4 

q(t) = q + —t- 1 - — r - — r + — — + ■■■ 

m I m 61 4! 5! 

-«°{ i -—+—+---)+^{ 9t -—+—+- 

= q cos 6t + -^-sm6t with 9 = \ — (12) 
mO V m 

which is fortunately the well known motion equation of a free harmonic oscillator. 

C. Systems with an infinity of freedom's degree 

The previous analysis was explicitly written for one particle (N = 1) but can be generalized to any N in a 
straightforward way. Let us study now the N — > oo limit, i.e., non-local statistical physics where variables are no 
more discrete. In this continuous limit, a dynamical system is no more described by a finite collection of {q s (A), p s (A)} 
but through a distribution function /(q, p, A) in the phase space. We deal with a non dissipative problem, hence, we 
suppose that the system evolve without collisions between its components. We will show that the structure of the 
equations and the solution of the problem of evolution of such systems are exactly the same than in the finite case. 

First, let us introduce a (non-relativistic) 4-dimensional formalism adapted to this study. We consider the 4— vectors 
Q = (q, A) and P = (p,fif), where g and A are the independents quantities defined in the previous section, which 
generalize energy and time for any kind of transformation. In the following we will consider functionals of/, considered 
as a function of {q, p, A}. In the space of functions of {q, p, A}, we denote the functional derivative b + . This notation 
is introduced to distinguish from the notation 6 for the functional derivative in the space of functions of {q, p} only, 
to be used later. We now define the meta-action S[f] through the relation j+y = g. In the special case where g does 
not depend on /, it can be written as 

S[f] := J f.g rf 3 q d 3 p d\, (13) 

but this linear formula does not apply, for instance, to the gravitational instability problem. We recall that we have 
assumed the generator g independent of A and we note that the meta-action S does not depend either on A. We also 
introduce another functional Q, now defined over the space of functions of {q, p} only (thus at given A), and such 
that |j = g. For the linear case where g does not depend on /, it can be written 
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<?[/] := f.g.d\d 3 p. (14) 



Note that, in general, Q depends on A. Let us define now a 4-dimensional Poisson bracket in the space of functions 
of {Q,P} by 

v *> y Pi[x ' y] = ^dQ.m-mm' (15) 

and a 4-dimensional Lie-bracket, relative to the distribution function /, in the space of functionals, as 

V X,y L 4 [X,y] (J) := j f.P.i—,^] (16) 

We are now able to show that the evolution equation of the system is equivalent to 

V/C L A [JC,S] (S) = 0. (17) 
Taking into account 8 + S[f]/8 + f = g, the latter relation is equivalent to 

V K, J f.P A [^j )9 ]d\d 3 v d\ = 0. (18) 



Following [7], an integration by parts (see appendix A) then gives 

6+IC 



K j 8 -^.PA[g,f]d\d^d\ 



0. (19) 



Since K, can be any functional, this last relation reduces to V^g, /] = 0. Developing 

and taking into account that g is explicitly independent of A, the equation is finally reduced to 

|{ = Ps[9Jl (21) 

This equation describes the evolution of the distribution function / of a non dissipative system. For the special case 
of motion, it is easy to check that it reduces to the usual Liouville equation. Making the mean field approximation, 
we recover the Vlasov dynamicalequation of collisionless systems. For continuous systems, observables are obtained 
by averaging physical quantities following the distribution function in phase space, thus observables may be seen as 
functionals K, acting on a positive and normalisable function / of {q, p, A}. Moreover for any such functional we have 

dfC[f] f 8fC df j3 j3 



Injecting (21) one has 

£ = /fftb,/K^P. (23) 

on inlo(VT"3fiAn Kir t~icit*+ lac 



Using g = |j, an integration by part (see appendix A) then gives 



dIC f r n r 6JC 8Q 1 3 



~x= ] fp * [ Tf>W** = Ms[fC ^ (24) 

where we have defined the new Lie brackets M3[JC,G](f) relative to the distribution function /. As for the discrete 
systems, M3[.,£J] is an A-independent operator, thus (24) can be solved using the Taylor expand of the exponential 

IC[f] = >C[fo] ~ A M 3 [G, JC] U) |a=a„ + y M 3 [G, M 3 [G, IC]] U) \x=x +■■■ (25) 
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For dynamical evolution (motion), the generator functional Q is the Hamiltonian 7i[f] of the system. Thus the 
time-evolution of any functional K, is given by 

£[/] = Wo] ~ t M 3 [H, JC] U) \t= to +^M 3 [H, M 3 [H, JC]] U) \t=t + • • • (26) 

This latter expansion gives the expression of any quantity (mean potential, entropy . . . ) of the system, at time t, 
as a function of initial conditions only. This may be applied to many problems, as far as the brackets are known. In 
the following, we apply it to the problem of gravitational instability. 

All these results can be resumed in a little box, expressing the symmetries between discrete and continuous non 
dissipative systems 



One degree of freedom 
Variables : q, p, A 

/c/q 
(ptt — g)dX 
dX 

Vfc(q,p) L 3 [k,s](q,p) = 

Evolution 



dk „ r , 
Solution : k(q, p) = e XPs[ ' 3] fc(q, p) | 



A = A„ 



Infinity of freedom's degrees 
Variable : /(q,p,A) 

M«,a-Ac„„„ : J ,,jA\,f vd X 

vicif] L,ir,s\„, = o 

Evolution 



fx = M ^ 

Solution : fC[f] = e XM ^ ^fC[f] | A=Ao 



III. GRAVITATIONAL INSTABILITY 



A. Dynamical equations 

Here we apply the previous results to the problem of Gravitational instability. This problem was introduced in this 
context by [4]. We will however not consider in this paper the standard formulation of the problem, using perturbed 
quantities, with respect to the Friedmann-Lemaitre, and a comoving formulation. This formulation would imply 
explicit time dependences in the evolution equations or the Hamiltonian, which would forbid the application of the 
previous rules. We rather calculate the evolution of the total (unperturbed plus perturbed) quantities. Since the 
evolution of the unperturbed ones are trivially known, under the form of the cosmological models, a simple difference 
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will provide the perturbations. We will consider the universe as filled of dust, i.e., pressureless matter only, and apply 
the Newtonian approximation as usual. Since gravitational instability becomes non linear much after the period of 
matter-radiation equivalence, i.e., when Universe is dominated by pressureless matter, this approximation applies 
well. 

The evolution of cosmic matter is described by the dynamical equation 
and the Poisson equation 

A J = 4ttG p, (28) 

where p is the matter density and <f> the gravitational potential. The velocity V = The time derivatives in 

these equations are to be taken following the motion, so that = J| + V.V q . Hereafter, primes will denote partial 
derivatives with respect to time J^-, i.e., derivative at a fixed (Eulerian) position q. The density evolves according to 
mass conservation p' = —div(p V). 

It is well known that it is possible to express the usual Friedmann-Lemaitre models in a (Newtonian) form obeying 
these equations. This will be our unperturbed solution : 

Pu = Po a~ 3 (29) 
V„=ifq= (a 1 1 a) q (30) 
K = 2ttG Pu q 2 A (31) 

where the scale factor a(t) obeys the Friedmann equations 

a 2 a" = -j/2 (32) 

and 

a' 2 = ja' 1 - k, (33) 

where we have defined j = 87r G po/Z. The scale factor a is normalized to its value at a time to, chosen as an origin 
and po denotes the value of the density at to ', k = —1, 0, or 1 is the curvature factor of space. 

B. The Vlasov— Poisson system 

Gravitational instability is usually described by hydrodynamic equations, to which the Poisson equation is added 
to express self-gravity. However, it is not always possible to find a Hamiltonian formulation of hydrodynamics. It is 
thus advantageous to work with the Vlasov-Poisson system, which is more general, and from which hydrodynamics 
can be deduced as a peculiar case. 

Thus we will consider equation (27) as the dynamical part of a Vlasov equation. This latter, which allows to follow 
the velocity distribution, is advantageous to consider for the Hamiltonian formulation. Hydrodynamics will be treated 
as a peculiar case of the Vlasov-Poisson formulation. Assuming particles with unit mass, the matter is described by 
a distribution function /(q,p,t) such that 

J (l 3 pp/(q,p,i), (34) 

and the average impulsion at point q 

p(q,t) V = J d 3 ppf(q,p,t) . (35) 

Note that the unperturbed distribution function takes the simple form / u (q, p,t) = po a(t) -3 6d(j> — H(t)q), where 
6d is the Dirac function. With these definitions, the Liouville equation takes the form 

^ +p ^_^^ = (36) 
dt 9q <iq dp 
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As we will check, the kinetic and potential "energies" associated to one particle are respectively T = p 2 /2 and </>(q). 
The gravitational potential at a point q is obtained by integrating (28) as 

<f,( qo ,t) = -G [ dT jf^'P'*) (37) 

J I q - q 1 

where dT = rf 3 pc? 3 p is the phase space volume element. 



C. Hamiltonian formulation of gravitational instability 

1. The Vlasov-Poisson system 

Considering gravitational instability as a motion with respect to time, its generator is the Hamiltonian g = T(p) + <f>, 
corresponding to the functional integral 

G = [ dTG/2 [ dT [ dr'lhi^lMiPlA (38) 

J J J I q - qo I 

so that the functional derivative |j = g. The application of the results of the previous sections involves the brackets 
^3[-4,0] (with respect to the function /). In order to calculate them, we need to use the following relations, 
straightforward to derive : 

¥UJ = ¥ =p (39) 



^ UfJ = ^ = ^ (40) 

where we wrote (/> q = . 

This allows to calculate the evolution of any functional A through the fundamental relation 

dA 

-£ = M 3 [A,g]. (41) 

Let us consider a functional A with its functional derivative jj- = A. When A does not depends on / (linear case), 
equation (41) implies that 

^ = j dTA^q^t) /(q,p,t), (42) 



where 



Al _ p dA dg_ dA dg__^dA dA 

9q '9p 9p '9q ' 9q ,q '9p' 



On the other hand, when A depends on / (non linear case), the situation is more complicated and, as we will see 

dt ' 



later, the functional derivative of A must appear in ^4- 



2. Hydrodynamics 



It is often sufficient - and this is the general approach - to consider gravitational instability from an hydrodynamic 
rather than Vlasov point of view. However there is no general Hamiltonian treatment for hydrodynamics (except in 
the case of irrotational flow). To recover hydrodynamics, we may start from the Vlasov description given above, and 
take the appropriate moments, neglecting the convenient terms, in the standard manner. In fact it is appropriate to 
remark that hydrodynamics for cold matter can be recovered by specifying the distribution function as /(q,p,t) = 
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p(q,t) 6d(p — ^(q,t))- Since the gas is considered to be cold and pressureless, all values of p at point q are equal to 
V. This allows to recover exactly the usual description of gravitational instability. 
The unperturbed distribution function for hydrodynamics is thus 

fu(<l,p,t) = p u (q,t) 6 D (p- V tl( q )t) ) (44) 

In the hydrodynamical case, the derivation rule derived above takes a simple form. When we are interested in the 
evolution of a functional quantity of the type 

A = j d 3 q p(q,t)A(q, V) , 

the application of the previous results leads directly to the derivation formula 

A' = J d 3 q p(q,t)A(q, V) , 

where A and A^ are related by (43). 



D. Time derivatives 



1. Time derivative of the distribution function 



It may be more interesting to evaluate the derivative of the gravitational potential at a point q . The expression 
(37) defines 4>(q ,t) as a functional of /, such that 



Ml) 



^(qo,0 



-G 



q - q 



Clearly, 



dp \6f 



(45) 



(46) 



dq. \6f) G dq VI q - q 



Applying the results of previous sections, we may now calculate the time derivative of <f> at point q 

^'(q0)0 = ~G J d 3 q d 3 p f(q,p,t) «i(q,p;q ), 



(47) 



(48) 



where 



«i(q,p;q ) 



G 



6<f> sg 
i] 'if 



p 



i 



dq V I q - qo 



Integration by parts, in the hydrodynamical approximation, leads to 

1 



b'(q Q ,t) = G / d 3 q 



div(p V). 



(49) 



q- q 

This equation could have been found directly by deriving (37) with respect to time and using the mass conservation 
= —div(p V) but here we illustrate how the method works to calculate any time-derivative. 
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2. Taylor development 



We have calculated the time-derivatives of two peculiar functionals, the distribution function itself and the potential. 
But we can also calculate, by applying iteratively the method, the derivatives at any order, of any functional A. This 
will allow, for instance, to calculate its Taylor development. For any functional A with functional derivative A, we 
have shown that the time derivative A 1 may be written as the integral of A^ related to A by (43). It results that 
the xi th order time derivative of A is of the form 

A {n) = [ dT f (50) 



where the A^ can be calculated from the A^ n This calculation involves the functional derivative 



(q, p) = A(")(q, p) + f dT' f ^"^''P') ^ (qj p) . (51) 



In this relation the last integral term takes into account the non linear case in which A^ depends on /. From the 
above calculations, it results that 



(52) 



which reduces to only 



^ (n+1) (q, P) = (P-^ - </v^) [^ (n) (q, P)L (53) 



in the linear case where A^ does not depend on /. 

Let us apply that to calculate some higher order derivatives of the potential. Let us first introduce the notation 
^(q — <lo) := |q-q | • Also, we will note the spatial derivatives any quantity A as 

A )8 j... fc = d ti d qj ...d tk A. (54) 

For instance, the J; = — jqZq"j3 are the components of a vector. With more indices, J 8 j...fc is a tensor. From 
c^q t -j = G f dT /(q, p, t) [—pi </,i(q — qo)L we can a PPly the previous formula to derivate again as 



nq„.*) 

with 



G I dT -^ q ' p '^ M2 ( q ' p;q o)' ( 55 ) 



«2 = - ( pj-; — ^>i~r~ j Pi J > 8 '( q ~ q °) = ~ Pi p i J > 8 j( q ~ q °) + ^ J > 8 '( q ~ q °)- ( 56 ) 

Reporting in the above formula, integrating by parts, and using the hydrodynamical approximation, we obtain 

1 ( d\p V, Vj] , d[p jf-] 



4" t) = -G rf d q 1 r L ; - Ji + • (57) 

(q °' ' J I q- q I V d<li dqj d qi J 

Subsequent derivation leads to higher order terms. Care must be taken however that the functional derivative of <f>" 
is not «2 which appears in the previous formula, since 1/2 is non linear and depends on the potential <f> which depends 
itself on /. Thus, to calculate from = <j>" , we must come back to original formulae (51). First we calculate 
the functional derivative 



"(q .«) 

6f 



( \ ( \ . f 7p/ fi ( ^ M 2(q',p';q ) ^ 

(q,p) = «2(q,p;q ) + J dT f I — — -l(q,p). 
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We continue by calculating the functional derivative of 1/2 as : 

jj ) (q, p) = J,,-(q - q ) I 6 j I (q)- 

After some calculations with no difficulty, the integral in the previous formula may be written 

(q'-qo)-(q-q') 



G j d\> „(q') j£ 



q I 3 |q-q' I 3 ' 



Finally, it results that 

^(q .*) 



0) c I d\ ( d*[ P V % Vj Vk ] | /[P d[ P V^] \ 

q - qo I \ <% <% d 1k dqi dqj dqj J 



+G [ d\> „ (q ') (q, - q ° ) -f q - q ?, 3 . (59) 

J I q - q 1 I q - q I 

Similar calculations would allow to calculate the next order terms and the Taylor development of the potential may 
be written as 

<f>(<lo,t) = <f>(<lo,t = 0) 

+G j rfr/ (q',p') (tu, + t 2 ^ + t 3 ^ + ...), (60) 

where all quantities are estimated at t = 0. 

There is no peculiar difficulty with the resulting formulae, at any order, excepted for their length, so that we will 
not write them explicitly. The potential at any point, and at any time (even after shell-crossing), may be calculated 
as a function of initial conditions only. But the expression is non local, involving the integration of initial quantities 
(density, velocity and potential) over the whole space. As expected, the expressions found are non linear, and it is easy 
to check that their degree increases with the degree in the time development. Thus the formulae here are comparable 
to the perturbative developments usually performed. Also it appears that each further degree in the development 
involves high order spatial derivatives of the initial physical quantities, integrated over space. This suggests that, the 
smoother the initial condition, the closer the solution remains from the linear one. 



IV. CONCLUSION 



We have presented a very general method to deal with Hamiltonian systems. Although the illustrating case is 
motion, i.e., evolution in time, this applies to any kind of evolution which can be described as a function of a time-like 
parameter, like invariances, instabilities, optics etc. . . For such systems, the evolution is described by a set of canonical 
variable (position and momentum for discrete systems), or by a distribution function in a continuous approach. 
Having introduced operators which express the algebraic structure of the problem, we have been able to present an 
evolution equation for any quantity depending on the canonical variables, or on the distribution function. We have 
also presented a formal solution of this equation, which allows to calculate any physical quantity as a function of 
"initial" conditions only, although under the form of a development with an infinite number of terms. 

These results apply well to the self-consistent problem of gravitational instability, where the distribution function 
evolves in time under its own gravitational interaction. This problem is treated here in the Newtonian approximation 
which applies in fact in the cosmological context, when fluctuations evolve with respect to a (dust) Friedmann- 
Lemaitre model. We present formulae which give an account of the evolution of any quantity depending on the 
distribution function (or on the mass density and velocity field in the hydrodynamic approximation). In particular we 
calculate the self-consistent evolution of the gravitational potential. We have expressed its value at any point, and 
at any time, under the form of a Taylor development whose terms of any order may be calculated as functions of the 
initial conditions only. This apply to the whole (unperturbed + perturbed) potential but it is easy to subtract the 
unperturbed part, which is known analytically, to derive the evolution of the perturbation. This will offer a practical 
way to calculate the evolution of the potential (or of any quantity like the density perturbations or the velocity fields), 
in the context of the gravitational instability theory for the formation of galaxies and large scale structures. 

In further work, we will compare the development introduced here with other development introduced in perturba- 
tive (Eulerian or Lagrangian) approaches. We will also treat in more details a case usually considered in gravitational 
instability scenarii, that with no initial velocity perturbations. In any case, this approach offers a new way to attack 
the question of gravitational instability. 
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APPENDIX : INTEGRATION BY PARTS AND SURFACE TERMS 



Let us consider 3 functions a;(q,p), j/(q,p) and ,z(p,q). We have 

/ Q3y Q3 Z Q3y Q3 Z 

\9q 3 9p 3 9p 3 9q 3 



j J xP 3 [y,z]d 3 q d 3 p = II . 



each term of the RHS of (Al) can be integrated by parts, 



— rf 3 q d 3 p 

aq 3 ap 3 



x y 



9p 3 



d 3 x d 3 z 



9p 3 9q 3 9p 3 9q s 



i 3 q rf 3 p 



(Al) 



(A2) 



for the first term and 



d 3 y d 3 z 
9p 3 9q 3 



f 3 q rf 3 p 



x y 



9q 3 



y 



( d 3 x d 3 z d 6 z 

+ x 

\9q 3 9p 3 9p 3 9q a 



i 3 q rf 3 p 



(A3) 



where Sq and Sp represents respectively the two 3-dimensional surfaces obtained when q and p goes to oo. Injecting 
(A2) and (A3) into (Al) one can get 



J I xV 3 [y,z]d 3 <i d 3 p = J I yV 3 [z, x]d 3 (\ d 3 p- 



dp 3 



9q 3 



(A4) 



In this last relation, the two last terms generally vanishes. Indeed, if 



d 3 z 
dp 3 



lim x y — — - = and lim x y -r-^- = 

— \3 :r-^. rin3 



d 3 z 
9q 3 



(A5) 



then these surface terms vanishes. For example if, x or y is the distribution function in the phase space of a Newtonian 
finite system, both this terms vanishes. For the unperturbed solution describing the cosmological models, it can also 
be checked that all surface terms also vanish, so that integration by parts can be applied as well. 
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